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Abstract 

We use the semiclassical method proposed in [1] to study scalar prod- 
ucts such as the overlap, Husimi functions and fidelity decay. Scars of 
classical periodic orbits arise naturally in this pertubative expansion. We 
also derive analytically a well known numerical result that fidelity has 
a quadratic decay for short times. We study the overlap in the chaotic 
regime and integrable for some simple systems. 

1 Introduction 

Since early times of quantum theory, some quantization difficulties of non in- 
tegrable systems were pointed by Einstein [2, 3]. Recently, due to the pioneer 
discoveries of classically chaotic systems, the subject has yielded many interest- 
ing and important results both from the point of view of numerical models and 
(not as many) analytical proofs [26, 31, 32]. 

One interesting discovery initiated by Bogomolny[4, 5] and Heller[6, 7] drew 
much attention. They showed that the Hamiltonians eigenfunctions of chaotic 
systems exhibit "scars" around the periodic orbit. The effect is more conspic- 
uous for integrable systems than for chaotic ones, given the wealth of periodic 
orbits in the latter [8]. The search for classical "imprints" than the celebrated 
phenomena of scars on eigenfunctions of quantum systems with classical analog 
has also gained a lot of attention. It was suggested [29] that a quantum spec- 
tra density and classical behavior are correlated . It was been demonstrated 
the existence of long-range correlation in quantum spectra and the existence of 
periodic orbits in the classical chaotic system [30] . Recently a relation between 
quantum phase transition and a classical instability points [9, 10, 11] was shown. 
All these cited works have a common characteristic, they show the existence of 
a close relationship between classical and quantum signatures. 

We begin by generalizing the semiclassical expansion [1] for n-dimensional 
system. The semiclassical expansion is builded such as the first order wave 
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function contains the classical dynamics for the system in question as completely 
as possible, in the sense that the dominant term is given only in terms of classical 
trajectories. All next to first order contributions contain essentially quantum 
effects and allow for a precise identification of departure from classical behavior 
in the quantum dynamics for short times. Secondly we use this expansion to 
extract classical ingredients known to be contained in several scalar products of 
states in the Hilbert space. We show in a very simple way how the phenomena of 
scars naturally arises from an adequate semiclassical analysis. On the other hand 
Husimi distributions and fidelity can be treated along the same lines and well 
known results from numerical studies can be given an analytical basis. The main 
advantage of the proposed method is due to his simplicity and his long domain 
of application, contrasting with the most actually used methods [26, 31, 32]. 



2 The Semiclassical Expansion 



Let us consider a classical one degree of freedom Hamiltonian of the form 

where p stands for the particle momentum and q for its position. We make a 
change of variables 

a — a* a + a* 
P= r— ; Q = 7= , (2) 

V mum y " 

where to = \Jkjm and k — 9 g i 9 ^ \ q =o- 

The Hamiltonian can then be rewritten as 

#ds = hhjaa + U(a*a) (3) 

with U(a*a) = V(q) - [k(a + a*)/ ^2muj/h} 2 . 
We can write H c \ s as a Taylor expansion, as 

ff c i s = frhj a* a + A m . n (a*) m a n , 

m.n 

where A\ t \ = 0. 

The classical equations of motions read 

m.n 

d „ 1 dH c \ s . i x ^ . m 1 , 

—a = - — — — = iuja + -> nA m . n {a ) a n \ (5) 
dt in da a ^ 



m.n 



We choose the quantum Hamiltonian in order to have (a\H q \a) = H c \ s , if \a) 
is a coherent field state. 



2 



We make our semiclassical expansion around a quantum operator H sc (a(t)). 
The difference, H — H sc (a(t)), will be considered as a perturbation .We choose 
the semiclassical hamiltonian, H sc (a(t)), that for a coherent initial state, all 
expectation values of point classical observables will be precisely reproduced. 

The semiclassical Hamiltonian which satisfies this condition is[l]: 

H sc (a(t)) =tkjtfa+J2 ™A m , m (a*)" 1 ' 1 a m -\a^a - a* a) + £ A m . n (a*) m a n 

m^O m.n 

+ Y, mA^ia*)™- 1 a n {tf - a*) + YnA m , n {a*) m a n -\a - a). 

m^n m.n 

(6) 

In the case of several degrees of freedom the semiclassical approximation is 
generalized as follows 

H sc (a(t)) = H cl + H£\a(t)) ®T 2 ■ ■ ■ ®Tjv 

+1 <8> flj? (a(t)) ® T 3 ■ ■ ■ (Sis + ■ ■ ■ 

Consider one degree of freedom we can immediately write the semiclassical 
evolution operator, just observing that 

U sc (t) \ao) = e im \a(t)) => U sc (a(t)) = e^V^^Ki))^)- 1 (7) 
where D stand for the well known displacement operator 

D{a)=e aa "- a ' a . 
Thus, for the N dimensional case we have 



U sc (a{t)) = Y[e* lt)i e a ' Wa i- a W a *e- a i a i +a i a i. (8) 

3 

N 
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The phase <p(a(t))j is given by 



cj>(a(t))j =~ij C[a,a*]dt' 



where C[a, a*] is the classical Lagrangian of the (independent) systems. In 
equation (8) we chose u)j(t) — 1 , since this can be done without violating any 



1 ojj refers to j-th system. 
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of the rules imposed to construct H sc {a(t)). It is just a matter of simplicity and 
can also be included, see [1]. In this case we are left with 



H(a(t)) = ^(a a(t)) - ^2 (tf _ a *(t)). (9) 

a stands for the on of the degrees of freedom in question. 

The generalization for SU(2) algebra or for any subspace where coherent 
states can be included, is immediate. The action of the semiclassical evolution 
operator on a coherent state can always be written as[12] 

uLWt)) |ft> = e ift(t) WAt) + ft - ^(o)) (io) 

where Pj(aj(t), ft) = (f> 3 (a{t)) - lm(aj(t)p*) + Im [(aj(0) - Pj)*a(t)] . In 
general case we have U 3 sc (a{t)) = £^ c (ai(t), ot2(t), ctz{t) 1 ajv(i))- 



3 Time Evolution 



Lets consider a two degrees of freedom system, whose complete Hamiltonian 
reads 

H = H 1 +H 2 + #12 (11) 

where #1(2) represent the autonomous dynamics of the degree of freedom 1 
(2) and Hi 2 is their interaction. The semiclassical Hamiltonian has the fol- 
lowing form H 8C (a,f3,t) = H] c (a,[3,t) + H 8C (a, fti), 2 and by definition we 

= 0. As discussed in section II we rewrite the 



have 



Hl c (a,(3,t),Hi(a,p,t) 



Hamiltonian (11) in the following form H = H sc (a,(3,t) + 5(a,(3,t), where 
S(a,(3,t) = H — H sc (a,(3,t) will be considered as a small perturbation. Using 
Schrodinger's equation we have 

ifiO. |M>) = H sc (a, ft t) I*) + 6(a, ft t) 

where we will always use |*(0)) = |a(0)) |/5(0)> , |a(0)) and |/3(0)> are coher- 
ent states. It is important to note that by construction we have 3 

(m\Ht)\d(a,MHt))\P(t))=0. 

Thus, after some straightforward [1] algebraic manipulations we get 



2 Hl c (a, f3,t) represents the semiclassical Hamiltonian in the subspace 1, and H sc (a, (3,t) 
represents the semiclassical Hamiltonian in the subspace 2. 

3 Since we have buildcd the semiclassical hamiltonian in such way that (^H^ = (^H sc ^. 
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I*) (t) = U sc (a,M { a 1 1 W * 

1 Jo Jo ^sit^A^dhdh 



l*> (0) 



(12) 



where A 8 (t) = Ul(a, /3,t)5{a, /3,t)U sc {a, /3,t). 



4 Scalar product of States 



Once working in a Hilbert space with a Hcrmitian Hamiltonian it is known 
that the scalar product of any two states which evolved with the same Hamilto- 
nian must remain constant. This "constancy" can be a test of our semiclassical 
approximation; moreover it teaches us something about when quantum cor- 
rections become dominant. The question we address is the following : let us 
consider two neighboring states in the sense that (f | $) w 1, wc know that 
| 4>) (t) = | $) (0). The important issue is: which of the ingredients of the 
quantum evolution are the ones that most affect this relation? Writing the states 
\*5f) and |<f>) in the form (12) it is a simple matter to obtain a semiclassical expan- 
sion for | $(*)) , for initial states |*(0)) = \a ,f3 ) and |$(0)) = \a' Q ,P' Q }. 
For short times we get 



where Uj c U' sc = U$ c {t,a,0)U sc (t,a',0'), e A\ = A s (ti,a,/3). Now intro- 
ducing the resolution of unity in terms of coherent states I = I\ <E> h = 



(*(t) | $(*)) = {p ,a \UlU' sc \a' ,l3' )-±- [ dti(foM&iUtcU'M,$>) 




J^J^\ti,(i) (C.^lweget 



<*(t)|$(t)) = <A )) ao|^ c 




+ 



The first term reads 



(/3 , ao| UlcU'sc K> P'o) = exp 



(\a t - a t '| 2 ) 



+ i Im (a t a' t *) + i{<f>(a'{t) - </){a(t))) 



2 



x exp 



(ift-fll 2 ) 
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We note immediately that this equation is relevant for the expansion for t 
where the conditions a t ~ a' t and fi t ~ (i' t hold. This t will be longer (almost 
always) in the case when the system is regular. 

Lets take a look on the time evolution of an operator a. We can compare it 
with his classical analogue, a c i. Suppose we can write the quantum expectation 
value as a function of his classical analogue, such as 

(a) q (t)=a cl {t) [1 + F q (t)] 

here, F q (t) is the quantum correction. Writing F q (t) 4 as 

- 1 d ™F q (t) 
q ( ) ~ ^ m! dt m 

m=2 

F q (t) then one defines the Erhenfest's time Time as 



The Erhenfest's time as it was described above, means the time where the 
quantum corrections are huge enough to become of the same order of the classical 
value. 

If we define a semiclassical time (r sc ) as the time during which we have 
| w (/?o, a \ UlcU'sc K, P' ) , then by definition [1] we have T E oc r sc , 

where Te being Erhenfest's time. Let us next consider systems with mixed 
dynamics, e. g. [14]. In this case we know that for the regular regions, neigh- 
boring trajectories keep close for a longer time than trajectories located in the 
ergodic region. The product ((3 ,ao\ U^ C U' SC \a' ,f3' ) will tend more rapidly to 
zero when we are in a region where there is chaos [14]. Figure (1) shows the 
semiclassical square modulus of the overlap between two neighbouring states, 
(ao \ U} C U' SC \ a' ) , for the driven conservative oscillator [15]. As we can see in 
this figure, the behavior of the semiclassical overlap is strongly dependent of 
the classical regime, as expected. Thus we can say that the validity of the 
semiclassical approximation is longer in the classical regular conditions. 

Consider two classical different initial conditions. Let D{t) be the distance 
of these trajectories in phase space. By definition, for t=T, D(T) >> DM =^> 
e -|£>(T)| <<; wnere dm = D(t sc ), i.e. it is the maximum distance in phase 
space that allows us to consider the first term as the most relevant. It means 
that there is an appreciable value for the product. Under these considerations, 
we can say that in the chaotic regime we have 5 r sc oc j ln(-^l), where A 6 
corresponds to the largest short time Lyapunov exponent. If we write Dq in 
terms of the canonically conjugate variables q,p then 



4 F q (t) This function in general depends on the initial state but is common to have a null 
first order. 

5 See Appendix. 

6 This exponent is calculated for short time series [13]. 
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(15) 
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We can thus say that Dq oc S c i, where S c i is the classical action of the system. 
This result is analogous to previously obtained ones [16, 17]. For regular regions 
we have a power law separation of neighboring trajectories, D(t) « D(0)t n , thus 



We can thus conclude that since the classical trajectories remain closer for 
longer times in the case of integrable systems, then we can say that they are 
more "robust" regarding quantum correlations, i.e., the following terms in the 
expansion. Chaotic systems on the contrary will very soon need quantum cor- 
rections for an adequate description. 



5 Correspondence Principle Aspects 

The phenomenon baptized as "scars" refers to hallmarks of the underlying clas- 
sical theory on its quantum counterpart, such as classical periodic orbits being 
very conspicuous in Husimi distributions. It was demonstrated that the exis- 
tence of scars can be shown by using semiclassical methods [6, 18], although 
the validity of these methods in the chaotic regime is not known, e.g. we know 
that WKB fails near caustic points[19]. In this section we will make use of the 
semiclassical expansion previously defined and show how it can shed light on 
the issue. 

5.1 Husimi's Quantum Phase Space Distribution 

The Q-function or Husimi's function is, see ref. [26, 27] , defined by: 



is a density operator, and |a) is the harmonic coherent state according 
to the definitions: 



H*(Q, P) = (a| | a) 



(16) 




(17) 



(18) 



From this definition, we are able to write the Husimi expression as: 
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H«(Q,P) = 



JJJ <f>(x,y,z)x(x,y,z)dxdydz 



(19) 



$(a;, y, z) is the studied system eigenfunction, and x( x i Vi z ) IS the harmonic 
coherent state in three dimension. x(x, y, z) can be written as: 



Xa(x,y,z) = Yl^i exp 



i=l 



(20) 



with 



6 4 W 



Mi 

Ax = Ay = Az = 



(21) 
(22) 



For the simplest case of the Harmonic Oscillator, using equation [16], the 
Husimi the Function for an eigenstate, n, can be written as: 



H„(a) = 



(27rfi)n! 



In terms of Q and P, we have 



H„(Q,P) = 



e -(Q 2 +P 2 )/2(Q2 +p2) « 



(47rfi,)n! 

5.2 Husimi Function for the Morse Potential 

The Morse potential is defined as : 

U(x) = D(e- 2ax - 2e- ax ), 

where we have defined 



(23) 



(24) 



(25) 



x = 



r -r 
ro 



(26) 



The r values are the equilibrium position of the center of mass, 
/i is the reduced mass of the two atoms. 
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The Hamiltonian that describes the center of mass can be written as: 



1 ,dr, 7 1 L 2 TT . . 
The time independent Schrodinger equation is: 



(27) 



"~ 1 '' 2 r+J^ + V(r)}<f>(r,6,(p)=E<S>(r,8,<p). 



2/i r dr 2 2\ir 2 
We can write the wavefunction as 



$(r,6,<p) = ^Y lm (6,<p), 



(28) 



(29) 



where Y/ TO is the spherical harmonics, so that is: 

Yr(e,v) = A lm Pr[co S (6)}e im «>. 

For L=0 case we find the eingcnvalues: 

E = -D + hm{(v+±)-±(v+±) 2 } 



(30) 



(31) 



and for the eigenfunctions: 



r (x + l) r(o) 



r(a + n) (C 2 e- a3: ) r 
r(<j + n) n! 



Where i/ = 0, 1, 2, 3, ... z/ < |(C 2 - 1), and 



(32) 



/3 2 = 


2^Er 2 

h 2 ' 


c = 


27 
a 


(2 = 


a 




f 1 


a = 


2 a 







Ai is fixed by normalization. 
Following the definition [21], we obtain the Husimi the Function as 

2 

H*(Q,P) 



/CO 
Jeo(a:)>r(a;)*(a;)ro(a; + l)dx 



(33) 
(34) 
(35) 



(36) 



where 7 eo (r) = -^ sinh( "° r) , 



0. 
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5.3 Semiclassical Husimi's Function 



The semiclassical expansion, as defined above, give us a time evolution of a 
quantum state as a pertubative expansion. A cigenstate has only a time de- 
pendent phase as its dynamics. The nearest semiclassical scenario we can build 
is to choose a coherent state with the same energy. The time dependence can 
be eliminated by a time integration, i.e. mean in time. This integration can 
be justified noting that as we are dealing with eigenstates we have not time 
precision. 

Under this considerations 7 we may write the semiclassical Husimi function 

as 



H ac (/3)« \<J3\a{t))\ . (37) 
The states \a(t)) and \/3) are coherent states of the harmonic oscillator. a(t) 

is defined as a(t) = anc i p = x + ip x , where x and p x are parameters 

of the Husimi Function. 

Q(t) and P(t) are the classical canonical conjugate pairs. 
Easily we can show that 



H sc (/3) w e \Q(t)+iP(t)-p\' ( 38 ) 
For the morse potential, with L=0, we obtain the classical trajectory [8]: 



1 , . D r . , 2 \ E \ a 2 . / D 2 + ED ^ 
a (t) = - M _ + Mt J-L^-W- [¥13 -]} E<0 (39) 



We also have p(t) = M^Ti and we can choose p(0) = and using E = E n 
into (38) to obtain the semiclassical Husimi function. 

In figure (2) we show the approximated Husimi for the Morse potential with 
the parameters of the Hi molecule, for n = 0. In figure (3) we have the exact 
result, figure (4) shows the semiclassical Husimi function for n=l and figure (5) 
the exact result, details about exact calculation can be found in ref. [8] . As we 
can observe in this figures, (2) and (3), the semiclassical Husimi function does 
not reproduce exactly the Husimi function, but it regards some similarities. 

Now consider the Harmonic potential, thus we have 

Q(t) = Q(0)cos(ojt) + ^^-sin(ujt) (40) 

and 

P(t) = P{<d)cos(ujt) - Q{0)ujsin(ujt) (41) 
We chose a coordinate system such as the Hamiltonian can be written as 



7 In case of classical mixed dynamics we must perform a mean considering all possible initial 
condition for the specific energy. 
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Substituting (40) and (41) into (38) we obtain a semiclassical Husimi Func- 
tion for the Harmonic oscillator with an energy E n = ftw{n + 5), without any 
lost of generality we can use Q(0) — (2E n )^ , and P(0) = 0. In figure (6) we 
show the approximated and exact Husimi for the Harmonic potential for n=5. 
Figure (7) shows the semiclassical and exact Husimi function for n=100. In 
the figures (3) to (9). In order to quantify the quality of the approximation we 
define the function AQ^ n (q) as 

AQ^ (q) = H(q,p=0) - H(q,p=0) sc , (43) 

where H is the exact Husimi function and H sc is the semiclassical Husimi func- 
tion. In Figure (8) we show AQ^^(q), for n=l and n=5. Figure (9) same graph 
for n=20 and n=100. Due to the spherical symmetry we have chosen p=0. X 
axis corresponds to position. As we increase the principal quantum number (n) 
we have AQ^ n (q) — > 0. In order to see the classical limit, let us define the 
function SQ^ n , witch is 

/oo 
\&Q^(q)\dq (44) 
-oo 

Suppose we have SQ^ « 0. It means that quantum description of the state tp, 
in the Husimi's representation, is almost contained in the semiclassical one. Of 
course it does not means that we have no quantum features, it only meas that 
Husimi is not a good observable for this situation [28]. Although that we can 
say that the quantum classical difference became smaller, as expected. Figure 
(10) shows SQ n for some eigenstates of the harmonic oscillator. For n=0 we 
have a null SQo, what was expected, since the fundamental harmonic oscillator 
eigenstate is a coherent state. Also we can say that the approximation works 
better as we increase the principal quantum number, as expected. From these 
figures we may conclude that the classical ingredient is very strong on the state 
formation of regular systems. 



6 Semiclassical Fidelity 

We now make a simple test of our approximation by applying it to a well known 
behavior of the fidelity. We know[20, 21, 22, 23, 24, 25] that for a linear pertur- 
bation , fidelity decays, for short times, as e~ at once one starts with Gaussian 

states. As we will see, this result is independent of the classical behavior of the 
system. 

Let us consider the product 

f{%V) = (HF\U {t)U{t)\9) (45) 

where Uq = e~% Ht and U = e~^^ H+v ^*, V is a perturbation, and it is 
given by V = eh(a + a* 1 ) where we have £ « 1. Fidelity (F) is defined as 
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F = \f(&, V)\ 2 . Let a(t) be the classical evolution of a(0) under the action of 
H , and (3(t) the classical evolution of [3(0) under the action of H = H a + V. 
We consider (3(0) = a(0). In this case for short times we have a(t) « (3(t) + iet. 
Our semiclassical expansion gives 

f(a,V) = (a(t)\\l3(t))e l ^ Wt)) -^ am 
1 /■« 

dtie 4 " 3 (a(ti)| tf[a(ti)] \a(h) - iet) 
' dtye^ (f3(h) + et\6[(3(h)} |/?(ii)> + -- ■ 



^ Jo 
^ Jo 



The term (a(i)| |/3(t)) gives the short time behavior for F and one gets 

e 2 t 2 



— iet 
2 



(46) 



(a(t)\ |/?(t)> = exp 

The above result is independent of the presence of chaos in the classical dy- 
namics. In Figure (11) we show the product \(ct(t) \ (3(t))\ sc . As we can easely 
observe, the chaotic initial condition decay faster. Although that, the short time 

dynamics has a gaussian decay, figure (12) we show \Jln\(a(t) \ f3(t))\ 2 sc for short 
time evolution, witch is a straight line. As the first term in the semiclassical ap- 
proximation is valid for longer times for intcgrablc systems than for chaotic ones, 
we conclude that the gaussian regime should be valid for longer times as it was 
already found numerically in several examples of the literature [20, 23, 24, 25]. 
Let us now look at the general initial state situation, f(a, V) can be written as 



<*| U (t)U(t) I*) = J ^ J *l m \a) (a\ U a (t)U(t) \(3) ((3\ |tt> (47) 

In the equation (47) we note that when a = (3, we have (a\ Uo(t)U(t) \(3) — 
f(a, V) as defined in (46), and it has the same characteristics discussed above. 
However, the product (*| |a) ((3\ (a\ Uo(t)U(t) \(3) is not expected to be very 
relevant when a is very different from (3, which is certainly true for short times. 
The whole analysis is weakened by the fact that the product ((3 \ may be 
important in chaotic and regular regions simultaneously. In this case it becomes 
very difficult to give a general estimate for the validity of (46). 



7 Conclusions 

As a general remark, we can say that, Overlap, Scar and short time Fidelity are 
strongly determined by the semiclassical dynamics, therefore we can say that the 
classical imprints are determinant. We also observe that overlap decay time has 
a very different dynamics from fidelity decay, although they are very similar in 
conception. In the particular harmonic case, we show that the first semiclassical 
term is able to reproduce the Husimi function with a increasing accuracy as 
we increase the principal quantum number n. We must remark that there is 
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no demonstration that would suggest an existence of a limit procedure witch 
turns quantum corrections less important in terms of the proposed semiclassical 
expansion. 
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Appendix 



A Semiclassical expansion and the Lyapunov Ex- 
ponent 

The semiclassical zeroth order is always a coherent state or a tensor product 
of coherent states,with the labels are described by classical dynamics. For a 
general state we have 

(tp\<j)} = {^\4>) sc + corrections. 

Consider that initially the states are \ip(0)) = JJ. |a 4 (0)) e \<j>(0)) = |A(0)> , 
where |oii(0)) e |/?i(0)) are coherent states. Then the zeroth semiclassical term 
is (ip(t)\(f>(t)} sc = n 4 (a*(i)|ft(i)) . We know [12]that the overlap of coherent 
states is 

|(a 4 ||/3 t )| 2 = exp(- |«i - A| 2 ) forh(3); (48) 



2 1 2J 

\oti - 0i\ 



for su(2). (49) 



(l + |ft| 2 )(l + k| 2 ) 

The equation (49) coincides with (48) if we make on = c^/v^J and taking 
the limit J — > oo. In this situation we can say that 

-|^-A| 2 = ln(|K|A)| 2 ) 

K-A| = v^anKailA)!" 1 ) 172 . (50) 
We can also say that \on — (3i\ 2 = (xi — yi) 2 + (p xi — p y i) 2 = Aq 2 + Ap 2 . 
The Lyaounov exponent is defined as 

A=lim Urn -In ^4 (51) 
*^»Ai(0)^0 t Ax(0) 

onde Ax(t) — \x\(t) — X2(t)\ , where Xi(t) is the clasical evolution for Xi(0) as 
initial condition. In the above limit, we get 

\{ at {t)\m)\ 2 = exp [- (Aq(0) 2 e 2 ^ + Ap(0) 2 e 2A « 4 )] 
As we have ^ A, = 0, the biggest Lyapunov exponent (A max ) is approximately 



Amax = lim lim — In 

t-ooft(0)- Q4 (0) It 



.InHli<«i(0)|A(0))r 
Observing this equation we may say that the Lyapunov exponent is related 
with the quantum nature of the system. As faster the quantum corrections are 
needed, i.e. how faster the product |(ct!j(£)|/3j(£))| — > , bigger is the Lyapunov 
exponent. This behavior has already been pointed by many others [17, 16]. 
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Figure 1: Squared modulus of the overlap between two neighbouring states 

2 

(oto\ Ul c U' ac \ot' ) , for the driven oscillator, x + x 3 — 0sin(u)t), with = 

1. Chaotic initial conditions (continuous line), lo = 1.88, (x\,xi) — 
(0, 0), (x 2 , x 2 ) = (0.002,0). Regular, same initial conditions (dotted line) for 
lo = 3.88. X axis corresponds to time. 
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Figure 2: semiclassical Husimi function for the Morse Potential, U(x) = 
D(e~ 2ax — 2e~ ax ), where x = r ~ r ° . We have used the experimental hydro- 
gen molecule values, a = 1.440 and D = 4.75(eV). The principal quantum 
number n=0 and the angular momentum L=0. The x axis is related to the 
position, r, as r = r o ^jgo > the momentum p is related to the y axis variable 
as p — mujr ^YY- . We have used r = 7.42 x 10~ n TO and u =8.3x 10 14 rad/s. 
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Figure 5: Husimi function for the Morse Potential with the same parameters of 
figure (4) 
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Figure 6: Full line shows a cross section of Husimi function for the harmonic 
potential for p=0. The principal quantum number n=5. The x axis is related 
to the position. Dotted line, correspondent semiclassical Husimi function. 




Figure 7: Full line shows a cross section of Husimi function for the harmonic 
potential for p=0. The principal quantum number n=100. The x axis is related 
to the position. Dotted line, correspondent semiclassical Husimi function. 
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Figure 8: Dotted line shows Qd function for the harmonic oscillator with the 
principal quantum number n=l, full line n=5. X axis corresponds to position. 
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Figure 9: Dotted line shows AQd function for the harmonic oscillator with 
the principal quantum number n=20, full line n=100. X axis corresponds to 
position. 




Figure 10: SQd function for the harmonic oscillator. X axis corresponds to the 
principal quantum number n. 
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Figure 11: Semiclassical fidelity \(a(t) \ 0{t))\ sc , for the driven oscillator, x + 
X s = jsin(u)t), with 7=1. Chaotic, u = 1.88, initial condition (dotted line) 
(xi,ii) = (0,0). Regular, lu = 3.88, same initial conditions (continuous line). 
Both curves with the same pertubation potential. X axis corresponds to time. 
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Figure 12: Scmiclassical overlap yln\(a(t) \ (3(t))\ 2 sc , for the driven oscillator, 
x + x a = jsin(cot), with 7=1. Chaotic, u) = 1.88, initial condition (dotted line) 
(x\,xi) — (0,0). Regular, ui = 3.88, same initial conditions (continuous line). 
Both curves with the same pertubation potential. 
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